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ABSTRACT 



BAT99-129 is a massive eclipsing binary system in 
the Large Magellanic Cloud (LMC) which consists of 
WN3(h) and 05 V components. A broad-band MACHO 
light curve of the system is studied in the present pa- 
per. A dense and extended atmosphere of the Wolf-Rayet 
(WR) star does not allow one to analyze the light curve in 
terms of standard parametric models of Wilson-Devinney 
type. Distributions of brightness and absorption across 
the WR star disk are restored using direct solution of in- 
tegral equations describing eclipses in the system. 

As a result, reliable estimates of the orbital inclination 
and component parameters have been obtained. The or- 
bital inclination is 78°, the orbital separation is 28.5 Rq, 
the radius of the O component is Rq = 7.1 Rq. The size 
of the WR component core which is opaque in optical 
continuum is i?vvR = 3.4 i?0. The brightness temper- 
ature in the center of the WR disk is Tbr = 45 000 K. 
Probable errors of the parameters are discussed. Evolu- 
tionary status of the system is discussed. 



1. INTRODUCTION 



Wolf-Rayet stars (WR) are considered to be helium rem- 
nants of evolved massive stars with initial masses of some 
dozens solar masses. A WR stage in the life of a mas- 
sive star does not last very long. Nevertheless, these 
objects play a key role in understanding evolution and 
internal structure of massive stars because of combina- 
tion of unique physical characteristics which make WR 
stars very valuable objects for comparison between the- 
ory and observations. One of the prominent features of 
WR stars is their strong stellar winds (mass loss rate may 
approach ^ IG^^'Mo/year), the origin of which is not 
fully clear yet. Dense winds make it very difficult to de- 
termine physical parameters of the central WR star cores. 
For example, two important parameters are the effective 
temperature of the WR star and the radius of its opaque 
core. This radius can be defined as a radius where the ra- 
dial Rosseland optical depth is equal to 2/3 1 1]. When the 
luminosity is known, the effective temperature of a star 
can be obtained using Stefan-Boltsman law. 



In principle, main characteristics of single WR stars (tem- 
perature, radius, and mass loss rate M) can be determined 
from comparison between their theoretical and observa- 
tional spectra. For dense winds, however, which are char- 
acteristic of early WN subclasses there is a correlation JH] 
between model parameters. The correlation allows only 
the combination of parameters L/ M^/^ to be determined 
where L — inRlaT^ff is the luminosity of the star, i?* is 
its radius. Toff is its effective temperature, a is the Stefan- 
Boltsman constant. 

Eclipsing binary systems containing a WR component 
provide a possibility to directly determine the radius of 
the opaque core and the brightness temperature of the WR 
by means of light curve analysis. A "normal" component 
in the system serves as a probe object. A strong semi- 
transparent wind from the WR makes parametric model- 
ing in binary systems (Russell-Merrill type |2], Wilson- 
Devinney type f3^, etc.) rather problematic. In this paper 
we use another approach, namely, direct solution of inte- 
gral equations describing light losses during eclipses. 

The second section of the paper contains the basic infor- 
mation on BAT99-129 and observational data which are 
used in our work. The third section describes the method 
of light curve solution. The forth section provides the re- 
sults of the light curve solution for BAT99-129. The fifth 
section presents a discussion of the results obtained and 
summarizes the conclusions of the work. 



INFORMATION ON THE SYSTEM. 
VATIONAL DATA 



OBSER- 



Eclipsing binary system BAT99-129 in the LMC is one 
of a few known extragalactic WR eclipsing binary sys- 
tems. Spectral observations of this star were recently an- 
alyzed in detail by |4]. Its spectral class was found to be 
WN3(h)H-05 V and the orbital period P = 2''. 7689 ± 
0''.0002. The orbit of the binary system seems to be cir- 
cular (zero eccentricity was accepted in the paper cited). 
The orbital separation is asini — 27.9 i?0. The ratio 
of the optical luminosities of the components has been 
determined using two spectrophotometric methods. The 
mean value of the ratio is L^f^/Lg"" = 0.34 ± 0.2. The 
radii and effective temperatures of the O and WR compo- 




Table 1. Mean Light Curve of BAT99-129. I - I is the 
light loss at a given orbital phase. 
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Figure 1. The light curve of BAT99-129 obtained in 
the MACHO experiment, individual measurements. The 
curve is reduced to relative intensities normalized to the 
maximum. The latter is determined as mean intensity be- 
tween the phases 0.3 and 0.35. 



nents were estimated approximately from their evolu- 
tionary status and corresponding model parameters. The 
effective temperature of the O component was accepted 
to be Tef / = 43 000 K while its radius is comprised be- 
tween 8.8 and 10.5 Rq. In reality, the O star may belong 
to 06 subclass as evidenced by some spectral features, 
thus having a smaller radius. The temperature of the WR 
component was estimated to be 71 000 K and the radius, 
which corresponds to the Rosseland optical depth of 20, 
i?* = 4.7 i?0. These values of radii as well as the sup- 
posed component masses allowed the authors to conclude 
that the inclination angle of the orbit may be around 60° 
or more. 
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BAT99-129 is an object in the sight of MACHO (Mas- 
sive Compact Halo Objectsfl In the course of the ex- 
periment a long-term photometric monitoring of a few 
million objects in the LMC and Milky Way is being car- 
ried out. The data are acqured in two color bands pro- 
dived by a dichroic filter The "blue" band is sensitive to 
emission between 4500-5900 Awhile the "red" one cov- 
ers the range 5900-7800 A. 

The data from more than seven years of observations 
for BAT99-129 system are available in both bands since 
1992. A total number of individual measurements 
reaches 877 in the blue band and 461 in the red one. Af- 
ter having been analyzed, the corresponding light curves 
do not show any significant differences in their shape. 
Hence, to analyze the light curve, the data in both bands 
have been combined and a joint light curve has been ana- 
lyzed instead. The individual measurements for BAT99- 
129 convoluted with the above mentioned period and the 
initial epoch [HJD] =2448843.8935 taken from 0] 
are shown in Fig.[T] 

The light curve in Fig.[T| is symmetrical with respect to 
the phase 0.5. For this reason, to carry out further anal- 
ysis, the right part of the light curve (the phases 0.5 - 
1 .0) has been reflected symmetrically relative to the phase 
0.5. After dropped-out points were removed and indi- 
vidual measurements were averaged, a mean light curve 

^The home page of the experiment is located at 
|http://wwwmacho.anu.edu.au| 



has been obtained in the phase interval 0.0-0.5. Inter- 
vals for the averaging were selected in such a manner 
that the mean light curve could describe light variability 
in the most optimal way. A mean-square error of a nor- 
mal point in the mean light curve lies between 0.002 and 
0.005. Normal points of the mean light curve together 
with their errors are given in Table[T] The table includes 
both non-rectified and rectified light curves (see below). 



3. METHOD FOR LIGHT CURVE SOLUTION 



To solve the light curve of BAT99-129 we used the 
method that was described in detail in [5] and [6]. Hence 
only its basic features are briefly described here and those 
who are interested in details are referred to the papers 
cited. 

Let us suppose that both components of an eclipsing bi- 
nary system are spherical and the brightness/opacity dis- 
tributions over their disks are axially symmetric. Then the 
light curve of the system containing a component with an 
extended atmosphere could be described by the following 
system of integral and algebraic equations: 



Ra. 

l-«i(6))= / A'i(e,A,i?o)/o/a(OdC (a) 



1-He)= / if2(C,A,iio)/c(0de {b) 



Rc 

LwR + Lo = h-J^RUl - |) + / /c(02vr^d^ = 1 



Here 1 — li,2{&) are observational light losses in the pri- 
mary and secondary minima respectively, is the rela- 
tive positional angle between the components connected 
with the distance A between the centers of the component 
disks by the formula A^ — cos^ i + sin^ i sin^ 6*, (the or- 
bital separation and the total luminosity of the system are 
suggested to be equal to unity), _Ro is the radius of the O 
star, Rc and Ra are the upper limits for the total radii of 
WR emitting and absorbing disks, /o is the brightness in 
the center of the O star disk, i is the orbital inclination, x 
is the limb darkening coefficient of the O star The func- 
tions Ki{£^, A, Rq), K2{(,, a, i?o) describe the shape of 
the overlapping disk area within the eclipses (their ana- 
lytical form is given in i?!], S)- lc{0 ^a(^) ™" 
known functions which describe brightness and absorp- 
tion distributions, respectively, over the WR disk. The 
latter can be represented in the form laiS.) = 1 — e"^'-^^ 
where t(^) is the optical depth along the line of sight in 
the extended WR atmosphere at the impact distance ^. 

Thus, the primary and secondary minima provide two in- 
tegral equations to be solved, their left-hand sides con- 
taining observational light losses and the right-hand ones 
involving the integrals of the unknown functions de- 
scribed above. As a matter of fact, the number of pa- 
rameters to be determined is not as large as it may 
seem. The brightness distribution over the disk of a main- 
sequence star can be satisfactorily described by a linear 
limb-darkening law. The problem thus involves only two 
respective parameters, the linear limb-darkening coeffi- 
cient, X, and the brightness in the center of the O star 
disk, Iq. Linear darkening coefficients for early main- 
sequence stars are well known from stellar atmosphere 
models; x = 0.3 for an 05 V star in the optical range. 
The central brightness Jo is determined from the third, 
algebraic, equation which connects the relative luminosi- 
ties of the components. Thus, the problem involves two 
unknown functions, la{0 ^i^d ^c(C)' together with two 
unknown geometrical parameters, the orbital inclination 
angle i, and the radius of the O star, Rq- It is important 
to note that one can only restore the brightness distribu- 
tion over the whole WR disk including its central parts, 
if the center of the WR disk is covered by the O star disk 
edge during the secondary eclipse. This condition can 
be expressed as an equation for the method applicability, 
cos(i) < Rq- We cannot know by default if this con- 
dition is fulfilled. There is, however, a simple and suffi- 
cient condition for the method applicability |10]: if the 
light loss in the lowest point of the secondary minimum 
(O star is in front) is equal to, or more than a half of the 



relative observational luminosity of the WR component, 
then clearly at least a half of the WR disk is invisible. If 
so, the center of the WR disk is certainly covered by the 
center of the O star 

The integral equations (la,b) are Fredholm equations of 
the 1st kind which are known to be ill-posed problems 
lUt]. Such problems could not be solved without some 
a-priori information concerning the behavior of the un- 
known function. The simplest example of such infor- 
mation is a parametric representation of the unknown 
function. For example, if a binary system consists of 
two main-sequence stars, the function IdO repre- 
sented by a linear limb-darkening law and /a(C) as being 
equal to 1 when ^ < i?wR and when ^ > i?wR- 
Then the problem could be reduced to a set of classical 
equations of the theory of eclipsing binary systems. 

For a WR star we can not use a linear limb-darkening law 
for /c(C)- Since the wind of the WR star is semi transpar- 
ent, a simple appproximation for {£,) given above is of 
no use either It was shown in [9] that an inverse prob- 
lem can be considered correct if its solution is searched 
on a compact set of functions. To solve such an inverse 
problem, standard methods can be used. Thus, to choose 
a type of a priori information acceptable in our case, one 
could ask a question: what minimal assumptions about 
the unknown functions can be done to make the problem 
correct in the classical sense? In |9], |8] it was shown 
that these assumptions can be very common in nature. In 
out case, the minimal requirements to the unknown Ic ) 
and la {^) are these: the functions must be restricted, non- 
negative and monotonically non-increasing. 

In practice these requirements tend to lead to solutions 
of low physical sense (for example H], @) when dis- 
tribution functions found are of a stepped form. Hence, 
in our work |6] dealing with light curve analysis of the 
eclipsing binary WN5 + 06 V V444 Cyg more strict re- 
quirements on the functions IdO and Ia{(.) were used. 
These functions were supposed to belong to a set of 
convexo-concave, non-increasing, monotonic and non- 
negative functions. A convex part of the functions cor- 
responds to the core of a WR star which contains most of 
stellar mass while a concave one describes the extended 
photosphere and the atmosphere of the WR star A po- 
sition of an inflection point is a free parameter of the 
problem. Such a priori information allows peculiarities 
of the WR star and its extended atmosphere to be taken 
into account in the least model-dependent way. The same 
requirements on the unknown functions are used in the 
current work. One additional condition is added, how- 
ever The radius of the opaque part of the WR core (where 
laiO — 10) is required to be equal to, or larger than the 
halfwidth of the function Ic{£,) at the half of its maxi- 
mum. In doing so, the widths of the functions IdO and 
la (C) are made consistent. Indeed, the central part of the 
WR disk should be absolutely opaque. If there were no 
such a restriction, a better description of the primary min- 
imum could be achieved for a function la {£,) hot having 
any flat part in the region of small ^. Typically, this hap- 
pens at large i. A "wide" flat-topped function Ia{£,) deep- 



ens the model primary minimum, and the only way for the 
minimization algorithm to reduce the divergence with the 
observational light curve is to decrease the width of the 
flat top in laiO- Because of that, the algorithm chooses 
as the "best" solution the one which does not have a flat 
top at all, i.e. the WR disk becomes semitransparent right 
from its center. The new restriction allows one to over- 
come this problem. 

Numerical solution of the problem (1) for a given pair of 
geometrical parameters i, Rq is done through discretiza- 
tion of the problem. The unknown continuous functions 
Jc(^) and laiO are represented as discrete functions de- 
fined in the knots of a rather dense grid over i^, the in- 
tegrals in (1) are replaced by sums, and the solution of 
the integral equations in (1) is found from residual mini- 
mization, the residual being a weighted sum of difference 
squares between the left and the right hands in equations 
(la,b). First the equation for the secondary minimum 
is solved, then Iq is determined from the normalization 
equation, and after that the equation is solved which de- 
scribes the primary minimum. 

Since the problem involves only two independent geo- 
metrical parameters (the orbital inclination and the O star 
radius) it is possible to look for the best solution by solv- 
ing the problem on a reasonable grid on these parameters. 
This method has the advantage of avoiding trapping in a 
local minimum. A normalized residual is defined as a 
weighted sum of squares of deviations between a model 
light curve and the observational one divided by the sum 
of weights. Since the problem is solved for two minima 
the resulting residual is a sum of two residuals obtained 
for the corresponding minima. 

The main disadvantage of the method is that it is impossi- 
ble to estimate how reliable is the resulting solution using 
standard statistical criteria of type- The matter is, the 
parameters of the problem are related in such a way that 
a real number of degrees of freedom could not be evalu- 
ated reliably. In theoretical works i^, |@] there is noth- 
ing but a mention that, while solving problem (1) on a 
compact class of functions, an approximate solution con- 
verges to the true one when observational errors decrease. 
In our case, we are dealing with a discrete light curve 
and a discrete representation for the unknown functions. 
The latters, generally speaking, not necessarily have to 
be smooth. Accuracy of the observational light curve as 
well as the number of data points it consists of are the two 
factors which define the quality of the approximate solu- 
tion. Indeed, if there was only one observational point 
in each eclipse the functions Ic{£,) ™d laiO could vary 
within huge limits still perfectly satisfying equations (1) 
and our a priori restrictions. For example, they can be 
mere constants. 

In our work (Tf] we carried out many numerical experi- 
ments solving artificially generated light curves by means 
of the above described method in order to find out the ef- 
fect of the mentioned factors on the solutions obtained. 
Here we briefly summarize the conclusions of this work 
in relation to an artificially created light curve that re- 



sembles the observed light curve of BAT99-129 (Fig.|2ll. 
The geometrical parameters of the system are determined 
very reliably, the deviations from actual values do not 
exceed 2-5 per cent. The radius of the opaque WR star 
core is also determined quite reliably (from the shape of 
the functions /c(0 ™d Ia{£,))- The value of /c(0) that 
is necessary to determine the brightness temperature of 
the WR star (see below) may have very large error. It 
is more likely to overestimate a real temperature than to 
underestimate it. A typical spread of the model /c(0) is 
—40% ~ +100% realtive to the "true" temperature. 

Using the function IdO obtained in solving (1) and as- 
signing the O star a value of effective temperature ac- 
cording to its spectral type, it is possible to determine the 
brightness temperature of the WR star in the central parts 
of its disk that is characteristic of the core temperature 



A = 
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(2) 



Note that the so obtained value of T;,(^, A) does not de- 
pend on interstellar absorption since the O star is used as 
a comparison star and the whole method of determination 
is differential. (For details see 0], 111, S). 



4. SOLUTION OF THE LIGHT CURVE OF 
BAT99-129 



The light loss in the center of the secondary minimum is 
1-/2(0) = 0.14. The relative observational luminosity of 
the WR component is — 0.25. Thus, the sufficient 
condition for applicability of the method is fulfilled: 1 — 

^2(0) > 5-/^wR- 

The light of BAT99- 1 29 between the eclipses exhibit reg- 
ular variability with the amplitude around 1 - 2 % (Fig.[T|). 
This means that the binary system components are not 
quite spherical. This may also testify about mutual heat- 
ing of the components. Our technique cannot take these 
factors into consideration. In the classical theory of 
eclipsing binary systems a rectified light curve is solved 
in such cases [2]. The standard rectification formulae (for 
example, for the model of rotating ellipsoids) cannot be 
applied to BAT99-129 (as one of the components is a WR 
star). We can derive, however, an empirical rectification 
formula approximating the light of the binary system out- 
side of eclipses: I — uq + ai ■ cos 9 + a2 ■ cos 29, where 
9 is a phase angle in radians |2]. A special question 
is revealing the phases which would signify the begin- 
ning and termination of the eclipses. Following recom- 
mendations from [2J, we have choosen the phase inter- 
vals 0.0-0.1 for the primary minimum and 0.4-0.5 for 
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Figure 2. Results of a test solution. An input light curve for given orbital parameters and functions /c(C)> -^a(C) 
calculated on the same orbital phases as the real observed light curve ofBAT99-129. Random deviations were then added 
to this "ideal" light curve distributed according to the normal law with root-mean- square deviation 0.003. In this way 100 
simulated light curves containg random errors have been obtained and solved. Black points stand for the true functions 
Ic{C) ond laiO- Solid lines represent the solutions for these functions for different realizations of the perturbed input 
light curve (20 solutions are shown to avoid encumbering the visualization). 



the secondary one. Approximation of the Hght between 
eclipses by the least-square method yields oq = 0.99311, 
= -0.01108, 02 = -0.00462. 

To better judge the reliability of our results we applied 
the above described technique to both non-rectified and 
rectified mean light curves of BAT99-129. The solutions 
are presented below. 



4.1. Non-rectified Liglit Curve 

The mean random error of a normal point in an aver- 
aged non-rectified light curve is 0.0030. The surface 
of residuals for the non-rectified light curve is shown in 
Fig.|3] Each pair of geometrical parameters Rq, i cor- 
responds to a specific relative luminosity of the compo- 
nents (see Fig.|4|i. The absolute minimum of the residuals 
{iji + il2)min = 0.0079 is achieved at Rq = 0.225, i — 
77° (marked by the cross in Fig|3] In this case, the rela- 
tive luminosity of the WR component is LwR = 0.377. 

As was mentioned in Section 2, the observational ratio 
of the binary component luminosities is L^'^^JLq^ — 
0.34 ± 0.20. This value corresponds to the relative WR 
component luminosity of WR L^f, = 0.25 ±0.11. The 
geometrical parameters corresponding to this luminosity 
ratio are shown in Fig. [3] as a solid broken line with 
short dashes. The minimum of the residual 771 + 1]2 along 
this line (marked by the triangle) is equal to 0.0080 (see 



also Fig. Is). The minimum corresponds to the parameters 
Ro = 0.25, i = 78°. 013. The residual in this point is 
practically identical to the residual in the absolute mini- 
mum and only slightly exceeds the doubled mean error in 
the observational light curve. Partly, this may be due to 
the fact that our model is not able to describe small regu- 
lar variability of BAT99-129 outside the eclipses. Hence, 
we have chosen, as a final solution, the parameters for the 
fixed WR component luminosity, iwR = 0.25: i = 78°, 
Ro — 0.25. Since the minimal residual exceeds the er- 
ror in the observational light curve, it is not possible to 
formally estimate (even roghly) uncertainties of these pa- 
rameters. The absolute orbital separation corresponding 
to this orbital inclination is 28.5 Rq. 

A model light curve and the functions lc{0^ ^a{C) for the 
chosen solution are shown in Fig.|6] The behavior of the 
function (^) indicates that the radius of the opaque WR 
core is 0. 13 of the orbital size (3.7 Rq). 

Using the value of /c(0) — 4.28, the parameters of the 
O star described in Section 2 as well as those obtained in 
the result of solving the problem (1), and the relative lu- 
minosity of the WR star, one can determine the brightness 
temperature in the center of the WR star disk by means 
of formula (2). This temperature specifies the core WR 
temperature. The wavelength in (2) is chosen to be equal 
to the mean wavelength of the MACHO spectral range, 
6150 A. The brightness temperature of the WR star core 
proves out to be 43 000 K. 
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Figure 3. The surface of the normalized residuals (bot- 
tom). The isolevels of this surface are shown on the top. 
The long-dashed line shows the boundary corresponding 
to the equation cos{i) = Rq. The short-dashed line cor- 
responds to a fixed relative WR luminosity iwR — 0.25. 
The cross marks the absolute minimum of the residuals. 
The triange marks the minimum along the line of the fixed 
WR luminosity. The isolevels which are the closest to the 
optimal solution have the values o/O.OlO, 0.015. 




Figure 4. The relative luminosity of the WR component 
versus geometrical parameters; the non-rectified light 
curve. 
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Figure 5. The residuals along the line of fixed luminosity 
Lwvi = 0.25 far the non-rectified light curve. 
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Figure 6. The solution of the non-rectified light curve for 
the optimal values of geometrical parameters. 

4.2. Rectified Light Curve 

The surface of residuals for the rectified light curve is 
shown in Fig.|7] The notations are the same as in Fig.|3] 
As in the case of the non-rectified light curve, the geo- 
metrical parameters of the solutions chosen according to 
the absolute residual minimum and according to the min- 
imum which corresponds to the fixed WR component lu- 
minosity, are slightly different. In both cases, the residu- 
als are practically identical, the absolute minimum being 
0.0053, and the residual at the fixed WR star luminos- 
ity, 0.0054 (see also Fig.[8]l. These values are somewhat 
lower than the doubled mean error of a data point in the 
mean light curve. As a final solution, we have chosen 
the one which corresponds to the fixed WR luminosity, 
iwR, = 0.25: i = 78°. 14, Rq = 0.25. 

In Section 2 we emphasised the difficulties to estimate a 
probable uncertainty of the solution. Very roughly, uncer- 
tainty in the geometrical parameters could be estimated 
by the isolevel which corresponds to the observational er- 
ror. In doing so, the observational error of i^f^ has also 
to be taken into account. The intervals 76°. <i< 79°. 5 
and 0.21 < Rq < 0.30 in the Rq, i plane define the area 
where the normalized residual does not exceed the obser- 
vational error and Lwr is constrained within the limits 
corresponding to the observational luminosity and its er- 
ror. 

The model light curve and the obtained functions IdO' 
la (0 for the chosen optimal solution are shown in Fig.|9] 
The opaque core radius of the WR star is 0.12 of the or- 
bital size (3.4 Rq). In this case, the WR brightness tem- 
perature is equal to 45 000 K (/^(O) = 4.52). 
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Figure 7. The surface of residuals for the rectified light 
curve. The designations are the same as in Fig.\3\ The 
isolevels which are the closest to the optimal solutions 
have the values of 0.006 and 0.010. 



One can estimate an uncertainty in the opaque core ra- 
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dius of the WR star and in the brightness temperature in 
the center of its disc. Let's find the i — Rq dependency for 
two fixed values of the WR star relative luminosity: 0.36 
and 0.14, which are the extreme values of the LwR un- 
certainty interval. Let's find the residual minimum along 
each i — Rq curve similarly to how this was done for 
LwK — 0.25. It turnes out that the two so obtained points 
on the Ro, i plane are located at the opposite sides of the 
error box given in the previous subsection. The solutions 
in these points yield the following values of the WR core 
radius and the WR temperature: 0.1 and 0.15 (2.8 and 
4.3 i?Q), 49 000K and 43 000 K, for Lwr = 0.14 and 
0.36 respectively. These limits can be accepted as rought 
estimates of the probable confidence intervals in the cor- 
responding parameters. 



5. SYSTEM PARAMETERS AND ITS EVOLU- 
TIONARY STATUS 



Based on our analysis of the MACHO BAT99-129 light 
curve the following conclusions concerning the parame- 
ters of the system and its components can be derived: 



Figure 8. The residuals along the line of the fixed lutni- 
nosity L^Y{. — Q.l^ for the rectified light curve. 
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1. The most plausible orbital inclination is 78°, a 
probable uncertainty interval being 76° — 79°. 5. 
Using the results from |4|]: asini = 27.9 i?©, 
MwR sin^ i = 14 ± 2Mq, Mq sin^ i = 23 ± 2Mq, 
we can estimate the absolute parameters of the sys- 
tem: a = 28.5 i?0, AfwR = 15 ± 2 Mq, Mq = 
24.6 ± 2 Mq. 

2. The most plausible O star radius is 0.25 of the dis- 
tance between the components, that is 7.12 Rq. A 
probable uncertainty interval for this radius is 0.21 - 
0.30, that is 6.0 — 8.6 i?© in absolute units. At an 
effective temperature of the O star 43 000 K its bolo- 
metric luminosity is Ig Lo /i© ~ 5.18. 

3. The most plausible opaque core radius of the WR 
star is 0.12 of the orbital separation (3.4 Rq). A 
probable error is 0.10 - 0.15(2.8 - 4.3 i?©). 

4. The brightness temperature of the WR core is 
45 000 K (its uncertainty is 43 000 - 49 000 K). Let 
us recall that this uncertainty is a rather formal one. 
Simulations show that a probable uncertainty can 
be up to some dozen per cent. It is more likely to 
overestimate the temperature than to underestimate 
it. The optimal temperature is in agreement with the 
results obtained from spectral modeling of isolated 
WR stars |1| and with our results obtained for the 
V444 Cyg [6] system. 



Figure 9. The rectified light curve solution for the optimal 
geometrical parameters. 



The function /a(C) can potentially be used to restore 
the spatial structure of the extended WR atmosphere |0]. 
This problem is also an ill-posed one. Such analysis 
was not carried out in present work because of great un- 
certainty of its results. To perform such analysis, high- 
accuracy narrow -band observations of BAT99-129 in op- 
tical and IR continuum are desirable. 



A small radius of the opaque WR star core (3.4 Rq), its 
mass equal to 15 AIq, as well as high brightness temper- 
ature of the core are indicative of an evolved massive hot 
star. The initially more massive component of a binary 
system can transform into a WR star in two ways (see, 
for example, |20]). It can lose the most part of its mass 
via stellar wind. The material leaves the system, the ac- 
cretion onto the companion being negligible. Both com- 
ponents evolve as if they were mere single stars. In the 
second scenario a more massive star fills its critical Roche 
lobe as it evolves and the so called initial mass exchange 
occurs in the system. The material flows from the more 
massive component to its companion. In wide binary sys- 
tems mass exchange is conservative and stellar material 
does not escape the system. In very close systems, a less 
massive component can also fill its Roche lobe after it 
accumulates enough material. Such a system becomes 
contact for a while. A common envelope forms around 
the system and some material can leave it taking away 
part of the system's angular momentum. 

The authors of |4] argue that the evolution of BAT99-129 
went through a contact stage. First, the WR wind mass 
loss rate is too small to explain the whole mass loss from 
the system on evolutionary timescale. Second, the period 
of the system (hence, its orbital separation) is very small. 
If mass exchange occured without the contact phase (a 
conservative case), the period and the orbital separation 
of the system would increase. The mass exchange sce- 
nario without contact phase implies that the initial period 
of the system should have been shorter than 2 days which 
is almost improbable. In contact mass exchange, how- 
ever, the orbital period and size are decreasing. Since the 
authors of |4] did not have precise information on the or- 
bital inclination of the system they suggested "typical" 
masses for the binary system components according to 
the spectral type - mass dependence. 

We have reliably estimated the orbital inclination of 
BAT99-129 which allows us to get some additional in- 
sight into the evolutionary scenario for the system. Let 
us suppose that initial conservative mass exchange took 
place in the system without a contact phase. In this case, 
the mass of a helium remnant (WR star) evolved from a 
massive star is related to the initial mass of its progeni- 
tor via the approximate relation ITtII : Mho — 0.1 • A/q'^. 
Using this formula leads to the initial mass of the WR star 
progenitor being equal to 36 AIq. Thus, the mass loss by 
the WR star during its evolution is about 21 A/q. Since 
mass exchange is supposed to be conservative the initial 
mass of the O component should have been only 4 Mq 
which is highly unlikely. These considerations show that 
the mass exchange could not be conservative, that is, 
some part of stellar material should have escaped the sys- 
tem. WR star formation according to the first scenario 
(mass loss due to stellar wind) appears to be unlikely, too. 
To transform into a WR star, the progenitor's mass should 
have been more than 25 Mq taking into account its rota- 
tion and metallicity Z — 0.008 lfl4ll . However, if a star 
had such initial mass it inevitably would pass through a 
supergiant stage before passing to the WR stage. In a 
close binary system such as BAT99-129, this would lead 



to contact mass exchange. A WR progenitor with greater 
mass (for example, > 50 Mq) which would imply the ab- 
sence of a supergiant stage in the course of evolution is 
less probable. In that case the system would inevitably 
lose a considerable part of material to form a circumstel- 
lar nebula. Nebula lines in the spectrum of BAT99-129, 
however, have not been found |4]. It looks liely that the 
initial mass of the progenitor amounted to no more than 
~ 40 Mq. a part of this mass has escaped the system at 
the stage of contact exchange. Nevertheless, this portion 
of mass is relatively small and has not been observed. To 
answer the question decisively one should obtain accu- 
rate spectra of the system with high spectral resolution 
and signal-to-noise ratio. 

While there appears to be some confidence in the param- 
eters and the evolutionary status of the WR component, 
the situation with the O component is far less certain. In 
Fig. [To] the position of the O star is shown on an evo- 
lutionary diagram for massive stars |15]. A dotted area 
stands for a probable error zone (see below). Evolution- 
ary tracks corresponding to the metallicity Z = 0.008 of 
the Large Magellanic Cloud (LMC) are used. Stellar ro- 
tation is not taken into consideration. The models with 
rotation are considered in other works (see, for example, 
||T6I] ). Unfortunately, evolutionary tracks with rotation 
have not been published for the LMC metallicity. In any 
case, a position of the Zero Age Main Sequence (ZAMS) 
on the evolutionary diagram does not depend strongly on 
rotation |16]. The star is located exactly on ZAMS. Its 
temperature and luminosity correspond to a ZAMS star 
with the mass of about 35 Mq, while the observational 
mass is only 25 Mq. Recall that the observational mass 
has been obtained using a rather reliable spectral estimate 
Mq sin^ i = 23 Mq |4] and our estimate for the orbital 
inclination. If the O star mass were 35il/0, the orbital 
inclination would have to be about 60 degrees which is 
absolutely incompatible with the results of our analysis. 

According to U5J, a ZAMS star with 25 Af© and Z = 
0.008 has the temperature 39 000 K and the luminosity 
IgLo/LQ = 4.87. Obviously, there exists an excess in 
the O star luminosity and its temperature. The mass and 
radius of the O component correspond perfectly to a the- 
oretical mass-radius dependence for ZAMS stars ( ITtII ). 
The position of the O component in BAT99-129 on this 
dependence is shown in Fig.[TT] The extra luminosity and 
temperature of the O component can be explained in three 
possible ways: 



1. Large errors in our T^ff and Lq. The error zone 
given in Fig. [TO] has been obtained in the following 
way. The error of T^ff is set to 1500 K (this er- 
ror defines the positions of the left and right borders 
of the error zone). This value is a typical intrinsic 
accuracy of temperature calibration which is set by 
comparison of atmosphere models with spectra of 
real stars. Unfortunately, the models used by dif- 
ferent authors yield systematic discrepancies which 
may considerably exceed the intrinsic accuracy of 
the individual models. As noted in lllSll . this may be 
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Figure 10. The evolutionary diagram for massive stars 
nlyil . Metallicity is Z = 0.008. To avoid cluttering, evo- 
lutionary tracks for the most massive stars are not shown 
completely. The bold point stands for the O component in 
the binary BAT99-129. The dashed line shows a probable 
error zone ( see the text). 
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Figure 11. The theoretical dependence of mass versus 
radius for ZAMS stars. The position of the O component 
in the binary BAT99-129 is shown. 



partly due to the fact that different models emphasis 
different spectral lines and spectral regions. In other 
words, contemporary models of stellar atmospheres 
are not quite adequate with real physics. As a result, 
"the scale of absolute effective O star temperatures 
is likely to be revised in future" |18]. Moreover, 
temperatures of different stars which belong to the 
same spectral subtype may differ significantly (even 
if a single atmospheric model is used [19]). Hence, 
the error zone shown in Fig.[TO]should be considered 
as minimal. On the other hand, the shortcomings 
of modern atmospheric models are likely to affect 
the ZAMS position on the evolutionary diagram as 
well. Thus, a relative position of the O star on this 
diagram would possibly not change too much when 
using another temperature calibration. 

The borders of the error zone defined by the error in 
the O star radius (the upper and the bottom borders) 
correspond to 6.0 and 8.6 Rq. Although it has been 
stressed above that we cannot provide truly sound 
statistical error estimates of the model parameters, it 
looks very unlikely that the radius of the O star is 
outside the above limits. 

Concluding, one can state that, while large errors in 
the parameters found cannot be formally ruled out 
as a reason for the "strange" position of the O star 
on the evolutionary diagram, this reason looks not 
very likely. 

2. The O component may belong rather to the 06 than 

05 subtype. This possibility is mentioned in H. 
The temperature and luminosity of an 06 compo- 
nent appear to be considerably lower than those of 
05. To coincide the position of the O component in 
BAT99-129 system in Fig[TO]with the initial point 
on the track for a 25 Af© star, its temperature should 
be To^cff = 39 000 K, and the radius Rq = 6.0 Rq. 
The temperature 39 000 K is typical for the Galactic 

06 stars and seems a bit low for the LMC ones (how- 
ever, see our notes above). The value of the radius, 
while being right on the bottom edge of a probable 
error box, seems to be somewhat too low. Finally, 
this opportunity can be verified when high-accuracy 
spectra of BAT99-129 are available which would al- 
low one to refine the spectral classification of the O 
star. 

3. In the process of the initial mass exchange, outer 
layers of the O component are enriched with the nu- 
clear burning products formed in the WR progenitor. 
This enforces mixing of stellar material [21]. The 
process is further supported by increased rotation 
velocity of the O component due to transfer of the 
angular momentum from the primary component. 
Spectral line widths of the O component in BAT99- 
129 reveal its possible asynchronous rotation fl. 
Thus, this component may be a chemically homoge- 
neous star enriched by the CNO-cycle products. The 
mass-radius dependence shown in Fig.[TT] is valid 
for any chemically homogeneous star [17], which 
may explain the position of the O component in this 
Figure. On the other hand, the anomalous chemical 



composition may be in charge of this component lu- 
minosity excess. Indeed, as was shown in |20], in 
massive MS stars in which the radiation pressure is 
dominant, the radius weakly depends on the chemi- 
cal composition. At the same time, other conditions 
being equal, the luminosity is proportional to the ra- 
tio of the mean molecular weight to 1 + X where 
X is the relative hydrogen abundance. Since the 
molecular weight increases and the value of X de- 
creases in the result of the initial mass exchange, the 
luminosity should somewhat rise compared to a star 
of the same mass and of the solar chemical compo- 
sition. 



6. CONCLUSION 



Wide-band observations of the eclipsing binary system 
BAT99-129 in optical continuum are interpreted in terms 
of an inverse problem solution. Specific a priori in- 
formation was used to constrain the unknown functions 
Ic (C) and la {£,) which represent brightness and absorp- 
tion distributions over the WR disk: they are convex- 
concave, monotonously non-increasing and non-negative 
functions. A convex part corresponds to the WR star core 
while a concave part describes the extended photosphere 
and atmosphere of the WR star The method used to solve 
the inverse ill-posed problem and to restore the unknown 
functions does not allow one to formally estimate the pa- 
rameter errors. To estimate the order of actual errors, and 
also the influence of the rectification procedure on the re- 
sults, both the original non-rectified and the rectified light 
curves were solved. As a result of such analysis the or- 
bital inclination of the binary system was estimated and 
the parameters of the components were determined along 
with their likely uncertainties. 

The analysis allows one to make some conclusions con- 
cerning the evolutionary status of the binary. The progen- 
itors of the binary components likely had masses around 
20^0 A/© and the system appears to have passed through 
a contact phase when it lost part of stellar material in the 
process of the initial mass exchange. The parameters of 
the WR star in the binary are typical for such objects. The 
O star shows some excess in its luminosity and temper- 
ature. The excess can be due to various reasons, most 
likely ones are uncertainty in spectral type determination 
and/or changes in the stellar chemical composition during 
the initial mass exchange. 

There are only a few eclipsing binary systems of the 
WRh-O kind where the parameters of the WR stars can be 
reliably determined. BAT99-129 is one of such systems. 
Its further observations in the optical and X-ray domains 
are highly desirable. 
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